Self-doping instability of the Wigner-Mott insulator 
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We present a theory describing the mechanism for the two-dimensional (2D) metal-insulator 
transition (MIT) in absence of disorder. A two-band Hubbard model is introduced, describing 
vacancy-interstitial pair excitations within the Wigner crystal. Kinetic energy gained by delocalizing 
such excitations is found to lead to an instability of the insulator to self-doping above a critical carrier 
concentration n = ric, mapping the problem to a density-driven Mott MIT. This mechanism provides 
a natural explanation of several puzzling experimental features, including the large effective mass 
enhancement, the large resistivity drop, and the large positive magneto-resistance on the metallic 
side of the transition. We also present a global phase diagram for the clean 2D electron gas as a 
function of n and parallel magnetic field B„, which agrees well with experimental findings in ultra 
clean samples. 



Significant experimental advances |1] over the last ten 
years have provided beautiful and convincing evidence 
for the existence of a sharp metal-insulator transition 
(MIT) in two-dimensional (2D) electron gases (2DEG). 
This progress has sparked much renewed interest in bet- 
ter understanding the basic physical mechanisms that 
drive the MIT, a fundamental physical question that has 
remained poorly understood for many years. 

One important issue relates to the stability of an in- 
teracting 2D metal with respect to disorder. While even 
weak disorder is known to destroy any 2D metal in ab- 
sence of interactions very recent work [s^l has provided 
strong theoretical evidence that electron-electron inter- 
actions may stabiHze a 2D metallic phase. This theory 
focuses on the most singular hydrodynamic corrections 
within the low temperature diffusive regime, and views 
disorder as the principal driving force that produces the 
insulating state. 

It should be emphasized, though, that the best evi- 
dence for a sharp MIT is found in the cleanest samples, 
where the diffusive regime is restricted to low densities 
and extremely low temperatures. Indeed, the experimen- 
tal data [l| demonstrating strong effective mass m* en- 
hancements have all been obtained in the ballistic regime, 
where diffusive processes are irrelevant. An important 
question then arises: How many of the key experimental 
features can be understood by deliberately disregarding 
disorder, and focusing on interaction effects alone as the 
main driving force for the transition? This question is 
the main subject of this paper, where we propose the 
quantum melting of a Wigner crystal as the fundamental 
mechanism for the MIT in a sufficiently clean 2DEG. 

It is well known that at the lowest carrier densities 
the 2D electrons form a triangular Wigner lattice. Here, 
each lattice site is occupied by a single spin 1/2 elec- 
tron, since strong on-site Coulomb repulsion prevents 
double occupation. The Wigner crystal can therefore be 
viewed as a (magnetic) Mott insulator, characterized by 
an abundance of low-lying spin excitations, but with an 
appreciable energy gap to charge excitations. As density 



increases, the gap to vacancy-interstitial pair formation 
decreases until the system undergoes a transition to a 
metallic state. The MIT from such a Mott insulator to 
a metal must, therefore, be fundamentally different from 
any Anderson-like transition, because the very physical 
nature of each insulating state is also completely differ- 
ent. If this idea is indeed correct - that the Wigner crystal 
melting is akin to a Mott MIT - then one may expect the 
critical behavior of the 2DEG to resemble that of other 
Mott systems. The canonical example for the Mott tran- 
sition in a continuum system is the Fermi Hquid to soHd 
transition in normal He^. Here, very recent experiments 
[4] on 2D monolayers have provided spectacular support 
for the Mott transition scenario. The effective mass m* 
was found Q to be strongly enhanced, while the g-factor 
remained essentially unrenormalized in the critical region 
- a behavior shockingly similar to that found in the most 
recent experiments [l| on the 2DEG! 

These arguments provide strong motivation to ap- 
proach the Wigner crystal melting as a Mott MIT, and 
develop an appropriate theory for the 2DEG. In the 
following, we describe the results of such an approach, 
demonstrating that the most striking experimental fea- 
tures of the 2D-MIT can all be simply understood within 
this framework. 

Charge transfer model. Our simplified description of a 
Wigner crystal is based on the idea that there exists a 
pronounced short range order in charge sector both on 
the metaUic and insulating sides of the MIT. This idea 
is strongly supported by quantum Monte-Carlo (QMC) 
work which shows that short range charge ordering 
changes little across the transition. It strongly suggests 
that, in a baUistic regime, a treatment in terms of an ef- 
fective lattice model should be applicable on both sides 
of the transition. The most important elementary exci- 
tations across the charge gap of a Wigner crystal corre- 
spond to vacancy- interstitial pair formation [g| • It should 
also be mentioned from the outset that our simplified de- 
scription does not treat collective modes explicitly, their 
effect being accounted only through renormaHzation of 
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effective microscopic parameters of the lattice. 

Sufficiently deep within the insulating phase, the elec- 
trons are tightly bound to lattice sites, and such excita- 
tions can be well described by an effective charge-transfer 
(e.g. two-band Hubbard ) model [t^ of the form 

■i(T ija 

+ E + 4/-) + E ufl^hflUi- (1) 

ia i 

Here, f\ f and , c are creation and annihilation opera- 
tors for site and interstitial electrons respectively, and U 
is the on-site repulsion preventing double occupation of 
lattice sites. In the tight-binding limit, the band struc- 
ture parameters t and V, and that of the charge-transfer 
gap Act = — £/, can be estimated by computing the 
appropriate wavefunction overlaps, leading to exponen- 
tially increasing bandwidth with density. The details of 
such band structure calculations will not affect any of our 
qualitative conclusions, and will not be reported here. 

The potentially most serious Hmitation of our lattice 
model is its phenomenological treatment (see below) of 
elastic deformation (e.g. collective charge excitations) 
of the Wigner lattice. These are expected to, at the 
least, effectively renormalize ^ the band structure pa- 
rameters, which should be quantitatively important for 
the physics of self-doping which we explore. More impor- 
tantly, one may question the very justification of using 
an effective lattice model, especially on the metaUic side 
of the transition. There is no particular reason, however, 
why the suppression of charge ordering associated with 
lattice formation must coincides with the closing of the 
Mott gap and the MIT. In absence of perfect nesting, the 
two transitions can occur separately, and the closing of 
the Mott gap may be expected to lead to an itinerant 
charge-ordered state. At any rate, the existence of short 
range order on the metallic side justifies, to a significant 
degree, the use of a lattice model in the inelastic regime. 

Mott transition via self-doping. In the Wigner-Mott in- 
sulator, the ground state the electron occupation is nat- 
urally one electron per unit cell. However, the lattice 
parameters in the 2DEG are self-consistently determined 
so that, as the density increases, it may become energet- 
ically favorable for the lattice spacing to slightly expand 
or contract, while keeping the charge density per unit 
area fixed (due to charge neutrality). If this happens, 
the resulting occupation per unit cell becomes 1 — 5, cor- 
responding to an effective self-doping of our Wigner-Mott 
insulator. Similar phenomena are believed to occur near 
the Mott transition in He'^ Q, and in the proposed for- 
mation of the He^ supersolid state Self-doping may 
change the precise nature of the MIT, and thus it needs 
to be carefully examined by properly incorporating the 
electrostatic considerations that are specific to a charged 
system of particles. 



Fermi-liquid condensation energy. Self-doping can be 
energetically favorable, since it leads to a kinetic energy 
gain of delocaHzed carriers which condense into a Fermi 
liquid. The price to pay is the cost of electrostatic energy 
to promote a carrier (electron or hole) across the charge- 
transfer gap. To assess the stability of the insulator to 
self-doping, one must calculate the doping dependence 
of the condensation energy of the incipient Fermi liq- 
uid state. This requires solving the appropriate Hubbard 
model - a general problem where no reliable or accepted 
theoretical approach is available at present. 

For our purposes a reliable treatment may be possible, 
and we seek inspiration from the closely related prob- 
lem of He^ monolayers Q. Here, the observed behav- 
ior can be quantitatively understood Q by the simplest 
Brinkmann-Rice (BR) theory of the Mott transition. 
This indicates that one approaches an insulator with lo- 
calized magnetic moments (hence a large m* as in any 
heavy fermion compound), where the inter-site spin cor- 
relations (measured by g*) in the Mott insulator can be 
ignored, as implied by the BR theory. Physically, this 
may be well justified for triangular lattices, where both 
the geometric frustration and the importance of ring ex- 
change processes (lol | conspire to render the spin correla- 
tions negligibly weak in the experimentally relevant en- 
ergy (temperature) range. 

To implement the BR approach for our problem, we 
follow the standard methods , where (for simplicity) 
we have taken U ^ oo. The free energy (per electron) of 
the self-doped system then takes the form: 



W[X,Z,fi,S] 
2T 



^^ln(l + exp(-(i?,fe-/i)/r)) 
A 



Ik 



1-5 



{z-^) + ^l (2) 



where T is temperature, Eik are renormalized band en- 
ergies, A is the Lagrange multiplier imposing the occu- 
pancy constraint, Z is the quasiparticle weight , and /i 
is the chemical potential. The free energy W[X, Z, fi,d] 
is stationary in the ground state: dW/da = 0, where 
a = X, Z, fi, S. 

The two bands of the model are coupled via the hy- 
bridization y/ZV. We further assume that the intersti- 
tial band density of states (DOS) is approximated by a 
regular function i/(e). The renormalized band energies 
then explicitly read : 



£^i,2(e) = ^{£f + £c + £' 



A) 



T^^Jiec-Sf + e-Xr+AZV^ (3) 



3 



The hybridization V and the density of states z^(e) have 
expHcit dependence on the lattice spacing and hence on 
the doping 5. The density depends on 5 as i^(e, S) — 
i^{e/j{S))/j{5), where ^{S = 0) = 1. The details of the 
dependence of V and 7 on i5 are not of qualitative impor- 
tance near the transition, as long as V{d) and j{S) are 
smooth functions of 6. Note that the van Hove singular- 
ities of a triangular lattice are sufHciently far from the 
Fermi energy at half filling. 

The choice of e/ and Sc, on the contrary, proves very 
important. We use effective electrostatic energy parame- 
ters to model complex energetics of the problem, arising 
from strong renormalization of bare parameters by elas- 
tic modes. For simplicity we assume a linear dependence 
of local potentials on the charge densities of the site and 
interstitial sublattices: 



(4) 



where = {/, c}, Uf = 1 — Z and Uc = Z — S, so 

the effective potentials depend on both S and Z. The 
prefactor — S represents overall rescaling of the 

Coulomb interaction with the change of the lattice spac- 
ing, due to charge neutrality. The coefficient Vij is an 
effective potential created on an (empty) site of the band 
i by the fully charged band j. We expect the polariza- 
tion to play a crucial role in the potential renormaliza- 
tion. When a hole is created the nearby electrons are 
attracted to the vacancy, partially screening it. When an 
electron is placed in the interstitial orbital, the nearby 
electrons are repelled again partially screening the 
charge fluctuation. The effect of the screening is al- 
ways directed toward decreasing the energy of the cor- 
responding particle-hole excitation. Due to the elastic 
softness of the Wigner lattice (e.g. shear phonons with 
the energy two orders of magnitude smaller, than the 
bare Coulomb energy) we expect a strong renormaliza- 
tion of the excitation energy, leading to an appreciable 
reduction in of the charge-transfer gap £c ~ £/ • There- 
fore we assume that (w// — Vcf)/vff — a \. The 
value of Vcc enters only the second order corrections in 
5, so we ignore it and simply set to zero. The stabil- 
ity requirement for the classical Wigner crystal restricts 
the value of Vfc- The electrostatic energy of the WC is 
E = (£/(! - Z) + £c{Z - 5))I{1 ~ 5). In the classical 
limit Z = (5 for (5 > and Z = for (5 < 0. By setting 
Vcf — {1 — a)vff we find that the WC is stable only for 

Solution of the model. The problem can be simplified 
in the critical regime (when Z ^ and (5^0) and 
analyzed analytically. This allows us to make certain 
general statement about the nature of the MIT. Away 
from the transition we resort to a numerical solution. 

a. Linear analysis and stabilization of the metallic 
phase via self-doping. In the following we show that the 
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Figure 1: (Color online) The phase diagram in the zero mag- 
netic field. The MIT takes route via self-doping (thick solid 
line), which always precipitates the transition at half-filling 
(thin solid line). The inset shows the quasiparticle weight, 
vanishing at ric. 



self-doping (SD) transition always precipitates the tran- 
sition taking place at half-filling (HF). We obtain a cri- 
terion for determining whether the SD is electron or hole 
like. 

We expand our equations around the insulating solu- 
tion {Z — 0, S = 0) to linear order in variations of the 
parameters A, Z, /i, S, assuming T — 0. At the transition 
point = + A and the free energy is purely classical: 
Wc — £f- From the saddle point equation dW/ dS — we 
find that in the SD transition W — fi = —dej/dS, hence 
Ac = dsf/dS. 

The results of the expansion in small Z can be con- 
veniently summarized in terms of the axillary functions 
(j)i{V,X) and 02(1^, A) defined as: 



dei>{e) 



1 



(f>2{V,X) = 1-2V^ 



Ef 



dev{e) 



Ec — £f + e — X 
(5) 

1 



iSa-£f+e- A)2 



(6) 



where Ep is the Fermi energy. One can check that the 
equation dW/dZ = for Z ^ (at SD or HF transi- 
tion) is simply (j)i{V, A) — 0. Differentiating this equation 
(where V — V{X) satisfies the equation) with respect to 
A we have: 



: 



d0i 

-dx=^' 



X+^dV 
V dX 



(7) 



It follows from the saddle point equations dW/ dX ~ 
and dW/dfi = that at the HF transition 02 (V", A) = 0. 
Therefore, according to the Eq.|[7l), dV/dX = at the 
HF transition. Direct inspection indicates that W has 



a maximum there, thus the SD transition always occurs 
before the HF transition. 

By considering dW/d\ — dW/d^i = near the SD 
transition one finds that Z — (5/(/)2(V^, Ac). Therefore, 
if (/)2 is positive in the SD transition, then the doping 
is hole-Hke, and if 4>2 is negative the doping is electron- 
hke. If 02 = the SD transition coincides with the MIT 
transition restricted to half filling. 

h. Numerical solution. We choose the parameters 
of the model that can best mimic the experimental re- 
sults. For that we set Vff — Ec, Vcf = (1 — a)Ec, 
Vfc = 0-lEc where a = 0.1 and Ec is the bare Coulomb 
energy. We use a parameter x = Dc/\Ec\ to mimic 
the Ts number, where Dc is the width of the intersti- 



tial band. The electron density goes as n 



in a 2D 



electron gas. We set V = Dc. The density of states 
in the interstitial band is constant, mimicking a two di- 
mensional dispersion. For these choice of the parame- 
ters we find (Eqs. (|5|6p ) that the SD transition occurs at 
xsD = 0.7408 {nsD = -549) and the HF transition occurs 
at XHF = 0.7751 {nsD = -601). 

In contrast to standard Mott transition, the half-filled 
insulator (heavy dashed line in Fig. 1) thus becomes 
unstable to electron-like self-doping (heavy full line in 
Fig. 1), before the half filled transition takes place 
(thin full fine in Fig. 1). The quasi-particle weight 
Z ^ 1/m* ^ {n ~ Uc) vanishes linearly (see inset of Fig. 
1) as the transition is approached from the metallic side, 
in agreement with experiments 

Transport, effect of magnetic field, and phase diagram. 
Different properties has been studied in detail for various 
Mott systems using recently developed DMFT method 
- a reliable tool even at low dimensions (unless the 
critical properties are specifically tied to the system's di- 
mensionality) . This approach can be regarded as a finite 
temperature generalization of the BR theory we utilized. 
Armed with this knowledge, one can directly list what 
is expected within the framework we consider: (1) Be- 
low the transition, transport takes place by activation 
p{T) exp (— A^(n)/T), with Ao{n) (x ric — n, just as 
seen in the experiments (2) On the metallic side, 
heavy quasiparticles exist only below a coherence tem- 
perature T*{n) ~ 1/m* ~ (n - ric), leading to a 
large resistivity drop [H] at T < T*{n); (3) A parallel 
magnetic field B*(n) ~ 1/m* {n — ric) is sufficient to 
produce full spin-polarization of the electrons, destroy- 
ing the coherent quasiparticles and causing large and 
positive magneto- resistance [ij; (4) Close to the tran- 
sition, at B» > B*{n) the resistivity saturates to a field- 
independent value p{T) — > poo{T), which assumes an 
activated form Poo{T) ^ exp (— A^(n)/r), where the 
gap Aoo(n-) remains finite in the high field limit. This 
behavior is specific to a charge-transfer (CT) model we 
consider, since the charge transfer gap Act remains fi- 
nite as — > oo, in contrast to the standard Mott gap 
Amou ^ U + gpsB,. (5) In the CT model the MIT 
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Figure 2: The metal-insulator phase diagram in the presence 
of a parallel magnetic field (in units of Be = ric/ f-iBgm). The 
dashed line represents a metal to band insulator transition. 
The electron spin becomes fully polarized at the magnetic 
field B* . The inset shows the change in resistivity from Fermi 
liquid to insulating behavior as the magnetic field exceeds B* . 



reduces to a band-crossing transition in the B„ — > oo 
limit, where Aoo(n) ~ (n^ — n) vanishes at > ric, 
and the system remains metaUic at higher densities. (6) 
The resulting phase diagram (Fig. 2) agrees well with the 
experimental one obtained for ultra-clean samples. 

Conclusions. We presented a theory for the 
interaction-driven MIT describing the clean 2DEG. 
Our approach focuses on vacancy-interstitial excitations 
within a Wigner-Mott insulator, naturally leading to a 
two-band (charge transfer) Hubbard model. As density 
increases, such excitations lead to an instability of the 
insulating state, and produce a self-doping driven Mott 
transition to a heavy electron metal. The general pre- 
dictions of this model seem to explain most puzzling fea- 
tures seen in the experiment, strongly suggesting that 
Coulomb interactions and not disorder provide the fun- 
damental driving force for the 2D-MIT. 

The most challenging task for future work is to extend 
the present approaches to explicitly include the dynamics 
of collective charge fiuctuations which are phenomenolog- 
ically treated in the considered lattice model. This goal 
should be facilitated by recent advances [l5| in theories 
for Coulomb gap phenomena, and would provide a more 
rigorous justification of the lattice model we introduced. 
Even more importantly, such a theory will be indispens- 
able to understand experiments at temperatures and 
densities where the Wigner lattice has already melted, 
but where strong short-range charge correlations per- 
sist. Such a regime is of appreciable importance and ex- 
tent whenever the reduced Coulomb interaction strength 
Ts 1, as found in many experiments on the 2DEG. 
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